Hadron structure effect in finite nuclei 
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The quark-meson couphng model, based on a mean-field description of non-overlapping nu- 
cleon bags bound by the self-consistent exchange of a, uj and p mesons, is reviewed. In 
particular, I present the changes in the hadron masses of the (non-strange) vector mesons 
and the nucleon in a nucleus. I also give a new, simple scaling relation among the changes of 
the hadron masses. 

§1. Introduction 

As the nuclear environment changes, hadron properties are nowadays expected to be mod- 
ified 11]. In particular, the variation of the light vector- meson mass is receiving a lot of 
attention, both theoretically and experimentally. Recent experiments from the HELIOS-3, 
CERES and NA50 collaborations at SPS/CERN energies have shown that there exists a 
large excess of the lepton pairs in central heavy-ion collisions [Q. On the other hand, even 
in low-energy nuclear physics, some attempts to measure the variation of hadron masses in 
nuclei are now underway [^. 

Recently we have developed the quark-meson coupling (QMC) model to treat the varia- 
tion of hadron properties in finite nuclei |^ . The QMC model may be viewed as an extension 
of QHD in which the nucleons still interact through the exchange of scalar and vector 
mesons. However, the mesons couple not to point-like nucleons but to confined quarks (in 
the bag). In studies of infinite nuclear matter it was found that the quark degrees of 
freedom in the nucleon give an acceptable value for the incompressibility once the coupling 
constants are chosen to reproduce the correct saturation energy and density for symmetric 
nuclear matter. This is a significant improvement on QHD at the same level of sophistication. 

In this report I briefly review the QMC model and show some numerical results for static 
properties and variation of the hadron masses in some nuclei |^] . I also give a new, simple 
scaling relation among the changes of the hadron masses in a nucleus. 

§2. The quark- meson coupling model 

Let us suppose that a free nucleon (at the origin) consists of three light (u and d) quarks 
under a (Lorentz scalar) confinement potential, Vc- Then, the Dirac equation for the quark 
field, il)q, is given by 

[i7 • a - mg - Vc(r)]V'g(r) = 0, (1) 
where niq is the bare quark mass. 



Next we consider how Eq.(|^) is modified when the nucleon is bound in static, uniformly 
distributed (iso-symmetric) nuclear matter. In the QMC model it is assumed that each quark 
feels scalar, V"/, and vector, V^', potentials, which are generated by the surrounding nucleons, 
as well as the confinement potential. Since the typical distance between two nucleons around 
normal nuclear density is surely larger than the typical size of the nucleon (its radius Rjy ~ 0.8 
fm), the interaction (except for the short-range part) between the nucleons should be colour 
singlet, e.g., a meson-exchange potential. Therefore, this assumption seems appropriate when 
the baryon density, pB, is not high. If we use the mean-field approximation for the meson 
fields, Eq.(|l|) may be rewritten as 

[ij-d- {m, - V!) -V,- = 0. (2) 

The potentials generated by the medium are constants because the matter distributes uni- 
formly. As the nucleon is static, the time-derivative operator can be replaced by the quark 
energy, —ieq. By analogy with the procedure applied to the nucleon in QHD, if we introduce 
the effective quark mass by m* = — V/, the Dirac equation, Eq.(|2|), can be rewritten in 
the same form as that in free space, with the mass m* and the energy tq — V^, instead of ruq 
and tq. In other words, the vector interaction has no effect on the nucleon structure except 
for an overall phase in the quark wave function, which gives a shift in the nucleon energy. 
This fact does not depend on how to choose the confinement potential, Vc- Then, the nucleon 
energy (at rest). En, in the medium is given by En = M'^iV^) + 2>V^, where the effective 
nucleon mass, M^, depends on only the scalar potential [0. 

Now we extend this idea to finite nuclei. The solution of the general problem of a composite, 
quantum particle moving in background scalar and vector fields that vary with position is 
extremely difficult. One has, however, a chance to solve the particular problem of interest to 
us, namely light quarks confined in a nucleon which is itself bound in a finite nucleus, only 
because the nucleon motion is relatively slow and the quarks highly relativistic [Q]. Thus the 
Born-Oppenheimer approximation is naturally suited to the problem. 

Our approach in Ref.|Q] was to start with a classical nucleon and to allow its internal 
structure to adjust to minimise the energy of three quarks in the ground-state of a system 
under constant scalar and vector fields, with values equal to those at the centre of the nucleon. 
In our case, the MIT bag model was used to describe the nucleon structure (see also Ref. [^). 
Of course, the major problem with the MIT bag (as with many other relativistic models of 
nucleon structure) is that it is difficult to boost. We therefore solve the bag equations in 
the instantaneous rest frame (IRF) of the nucleon - using a standard Lorentz transformation 
to find the energy and momentum of the classical nucleon bag in the nuclear rest frame. 
Having solved the problem using the meson fields at the centre of the nucleon, one can use 
perturbation theory to correct for the variation of the scalar and vector fields across the 
nucleon bag. In first order perturbation theory only the spatial components of the vector 
potential give a non- vanishing contribution. This extra term is a correction to the spin-orbit 
force lH]. 

As shown in Refs.|^, H], the basic result in the QMC model is that, in the scalar and vector 
fields, the nucleon behaves essentially as a point-like particle with an effective mass M]^, 
which depends on the position through only the scalar field, moving in the vector potential. 

Let us consider that the scalar and vector potentials in Eq.(^) are mediated by the a 
and uj mesons, and introduce their mean-field values, which now depend on position r, by 
V^{r) = g%(j{r) and V^{r) = g^uj{f), respectively, where g$ (g-^) is the coupling constant 
of the quark-cr (a;) meson. Furthermore, we shall add the isovector, vector meson, p, and 



the Coulomb field, A(f), to describe finite nuclei realistically |Q, ^. Then, the effective 
lagrangian density for finite nuclei, involving the quark degrees of freedom in the nucleon and 
the (structureless) meson fields in MFA, would be given by 



C-QMC-I = 



- ^[(Va)2 + rn2a2] + i[(Vc.)Vm2^2j^i[(V6)2 + m262] + i(vA)2, (3) 

where ^ip and b are respectively the nucleon and the p fields, mg-, and nip are respectively 
the (constant) masses of the u, uj and p mesons, and gp are respectively the w-N and 
p-N coupling constants, which are given by 5^ = 8(7^ and gp = g^, (where g^, is the quark-/? 
coupling constant). We call this model the QMC-I model. If we define the field-dependent 
(T-N coupling constant, 5o-(cr), by 

M*Mr)) ^Mn- g.{a{r))a{r), (4) 

where M^r is the free nucleon mass, it is easy to compare with QHD. The difference between 
QMC-I and QHD lies only in the coupling constant g^, which depends on the scalar field 
in QMC-I while it is constant in QHD. However, this difference leads to a lot of favorable 
results, notably the nuclear incompressibility [Q]. 

Here we consider the nucleon mass in matter further. The nucleon mass is a function of 
the scalar field. Because the scalar field is small at low density it can be expanded in terms 
of a as 

M'„ = M„ + (^) , + - (^ j + . . . . (5) 

Since the interaction Hamiltonian between the nucleon and the a field at the quark level 
is given by Hint = — / df ipqaipq, the derivative of with respect to a is given by 
—Sg'^ J df tpijipg(= —3g^SN{o')), where we define the quark-scalar density in the nucleon, 
Si\i{a), which is itself a function of the scalar field. Because of a negative value of the 
derivative the nucleon mass decreases in matter at low density. 

Furthermore, we define the scalar-density ratio, SN{cr)/ S]y{0), to be Cjyio') and the cr-N 
coupling constant at cr = to be g^ (i.e., CNicr) = 5Ar(o")/5'Ar(0) and ga = 5o-(o" = 0) = 
3g^SN{0)). Using these quantities, the nucleon mass is rewritten as 

M*^ = Mn- ga(T - ^gAma^ + ■■■. (6) 

In general, C^r is a decreasing function because the quark in matter is more relativistic than 
in free space. Thus, C^(0) takes a negative value. If the nucleon were structureless Cn 
would not depend on the scalar field, that is, Cn would be constant (Cat = 1). Therefore, 
only the first two terms in the right hand side of Eq.(|^ remain, which is exactly the same 
as the equation for the effective nucleon mass in QHD (see also Ref. ^^). By taking the 
heavy-quark-mass limit in QMC we can reproduce the QHD results ^, ^. 

We have considered the effect of nucleon structure. It is however true that the mesons are 
also built of quarks and anti-quarks, and that they may change their properties in matter. 
To incorporate the effect of meson structure in the QMC model, we suppose that the vector 
mesons are again described by a relativistic quark model with common scalar and vector 
mean- fields like the nucleon (see Eq.(Q)). Then, again the effective vector-meson mass in 
matter, m* {v = co 01 p), depends on only the scalar mean-field. 



However, for the scalar (cr) meson it may not be easy to describe it by a simple quark model 
(like a bag) because it couples strongly to the pseudoscalar (27r) channel, which requires a 
direct treatment of chiral symmetry in medium [0]. Since, according to the Nambu~Jona- 
Lasinio model or the Walecka model, one might expect the u-meson mass in medium, m*, 
to be less than the free one, we shall here parametrize it using a quadratic function of the 
scalar field (by hand) [^: 



^) = l-aM„a)+b^{g„ay, (7) 

with ga-(7 in MeV, and we introduce two parameters, Og- (in MeV~^) and 60- (in MeV~^). We 
will determine these parameters later. 

Using these effective meson masses, we can find a new lagrangian density for finite nuclei, 
which involves the structure effects of both the nucleon and mesons in MFA: 



C-QMC-II = i> 



^ (8) 



We call this model QMC-II. Note that all the masses depend on the scalar mean-field. 



§3. Numerical results 



In this section, we will show our numerical results using the QMC-II model. (For QMC-I, 
see Ref. [|.) 

For infinite nuclear matter we take the Fermi momenta for protons and neutrons to be kp^ 
(i = p or n). This is defined by pi = kp./{3TT^), where pi is the density of protons or neutrons, 
and the total baryon density, /9_b, is then given by pp + pn- (Note that the mean- field values 
for the mesons are constant.) 

From the lagrangian density Eq.(P), we can easily find the total energy per nucleon, Etot/A, 
and the values of the uj and p fields (which are respectively given by baryon number conser- 
vation and the difference in proton and neutron densities, Ps = Pp — Pn)- On the other hand, 
the scalar mean-field is given by a self-consistency condition (SCC): = 0. 

Now we need a model for the structure of the hadrons involved. We use the MIT bag 
model in static, spherical cavity approximation |jl2|. In the model, the bag constant B 
and the parameter zat in the familiar form of the MIT bag model lagrangian are fixed to 
reproduce the free nucleon mass (M^v = 939 MeV) under the condition that the hadron 
mass be stationary under variation of the free bag radius {Rn in the case of the nucleon). 
Furthermore, to fit the free vector-meson masses, = 783 MeV and mp = 770 MeV, we 
introduce new z-parameters for them, and Zp. In the following we choose Rn = 0.8 fm 
and the free quark mass niq = 5 MeV. We then find that B^^'^ = 170.0 MeV, zn = 3.295, 
Zo; = 1.907 and Zp = 1.857. (Variations of the quark mass and Rn only lead to numerically 
small changes in the calculated results p.) 

Next we must choose the two parameters in the parametrization for the a-meson mass in 
matter (see Eq.(^)). Here we consider three parameter sets: (A) Oo- = 3.0 x 10~^ (MeV~"^) 
and ba = 100 x 10"^ (MeV^^), (B) a„ = 5.0 x lO"^ (MeV^^) and 6„ = 50 x 10"^ (MeV'^), 



(C) = 7.5 X 10"^ (MeV-i) and b„ = 100 x 10"*^ (MeV-^). The parameter sets A, B and 
C give about 2%, 7% and 10% decreases of the a mass at saturation density, respectively [^. 

Now we are in a position to determine the couphng constants, and are fixed to fit 
the binding energy (—15.7 MeV) at the saturation density (po = 0.15 fm~^) for symmetric 
nuclear matter. Furthermore, the Qp is used to reproduce the bulk symmetry energy, 35 MeV. 
We take nia = 550 MeV. The coupling constants and some calculated properties for matter 
in QMC-II are listed in Table ||. The last three columns show the relative changes (from their 
values at zero density) of the nucleon-bag radius [5R*^ / Rn), the lowest eigenvalue {5x*^/xn) 
and the root-mean-square radius (rms radius) of the nucleon calculated using the quark wave 
function {6r*/rq) at saturation density. We note that the nuclear incompressibility is higher 



Table 1: Coupling constants and calculated properties for symmetric nuclear matter at normal 
nuclear density {niq = 5 MeV, Rn = 0.8 fm and nia = 550 MeV). The effective nucleon mass, 
M^, and the nuclear incompressibility, K, are quoted in MeV. The bottom row is for QHD. 



type 




5i/4vr 


5^/4vr 




K 


5R*j^/ Rn 




6r*/rg 


A 


3.84 


2.70 


5.54 


801 


325 


-0.01 


-0.11 


0.02 


B 


3.94 


3.17 


5.27 


781 


382 


-0.01 


-0.13 


0.02 


C 


3.84 


3.31 


5.18 


775 


433 


-0.02 


-0.14 


0.02 


QHD 


7.29 


10.8 


2.93 


522 


540 









than that in QMC-I {K ~ 200 - 300 MeV) [|, However, it is still much lower than in 
QHD. As in QMC-I, the bag radius of the nucleon shrinks a little, while its rms radius swells 
a little. On the other hand, because of the a field, the eigenvalue is reduced more than 10% 
(at po) from that in free space. 

In the present model it is possible to calculate masses of other hadrons like hyperons in 
medium (see also Ref. [^].) In that case, we assume that for the hyperons themselves we 
again use the MIT bag model, and that the strange quark in the hyperon does not directly 
couple to the scalar field in MFA, as one would expect if the a meson represented a two- 
pion-exchange potential. It is also assumed that the addition of a single hyperon to nuclear 
matter of density pB does not alter the values of the scalar and vector mean- fields, namely, 
we take the local-density approximation to the hyperons. 

In general, as in the nucleon case (see Eq.(|5[)), we can expand the effective hadron mass in 
nuclear matter at low pB as 




M,-^g.cT-^5.q(0)a^ (9) 

where j stands for A^, co, p, A, E, H, etc., no is the number of non-strange quarks in the hadron 
j and Cj{a){= Sj{a)/ Sj{0)) is the scalar-density ratio defined by the quark-scalar density, 
Sj, in j. Since the strength of the scalar mean-field, g^^a, at low pB can be approximated by 
a linear function of the density |5[ as g^a ~ 200 (MeV) {pb/ Po)} we find the hadron mass at 
low pb'- 




where cn,v,a,s,s = 0.21,0.17,0.12,0.11,0.05, respectively . 

Furthermore, using the MIT bag model we find that the scalar-density ratio Cj can be 
well approximated by a linear function of the scalar field |||]: Cj{a) = 1 — aj x (^o-c), where 
aj is the slope parameter for the hadron j, and that the dependence of Uj on the hadrons is 
quite weak (it ranges around 8.6 ~ 9.5 x 10"^ (MeV"^)). Therefore, if we ignore its weak 
dependence on the hadrons, the effective hadron mass can be rewritten in a quite simple form 
(from Eq.@): 

no ^ , r. a . 

Y 



1 - i^{9a<y) 



(11) 



where a ~ 9.0 x 10~^ (MeV"^). This mass formula can reproduce the hadron masses in 
matter quite well over a wide range of pB (up to ~ 'ipo). Therefore, once one knows the 
behaviour of g^a at finite density, one can easily calculate the effective hadron mass using 
Eq.([Tl|). A simple parametrization of gfjG is given in Ref. fl^ . 

Since the scalar field is common to all hadrons, Eq.([lTl) leads to a new, simple scaling 
relationship among the hadron masses |5[: 

5mt\ ( 8Ml\ ( 5M^\ 2 ^ ( 5Mt\ 1 

and — - o' (12) 



J \6M^ J \6M^ J 3 \SM^ J 3' 

where 6M* = Mj - M*. The factors, § and i, in Eq.(0) come from the ratio of the number 
of non-strange quarks in j to that in the nucleon. It would be very interesting to see whether 
this scaling relationship is correct in forthcoming experiments. 
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Figure 1: Changes of the nucleon, a and uj meson masses in Ca. The nuclear baryon density 
is also illustrated (solid curve). The right (left) scale is for the effective mass (the baryon 
density). The parameter set B is used. 



For finite nuclei, we will show our results of some finite, closed shell nuclei. We have solved 
self-consistently a set of coupled non-linear differential equations given by the lagrangian 
density Eq.(^. In Fig. |l] we present our self-consistent, full results of the changes of the 
nucleon, a and lo meson masses in '*°Ca. Table ^ gives a summary of the calculated binding 



energy per nucleon (E/A), rms charge radii and the difference between nuclear rms radii for 
neutrons and protons (r„ — Vp) for several closed-shell nuclei ^. While there are still some 



Table 2: Binding energy per nucleon, —E/A (in MeV), rms charge radius Vch (in fm) and the 
difference between r„ and (in fm) for type B, niq = 5 MeV and Rb = 0.8 fm. I and II 
denote, respectively, QMC-I and QMC-II. (* fit) 





-E/A 


rch 


^p 


model 
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II 


Expt. 
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II 


Expt. 
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II 


Expt. 


ItiQ 


5.84 


5.11 


7.98 


2.79 


2.77 


2.73 


-0.03 


-0.03 


0.0 


40Ca 


7.36 


6.54 


8.45 


3.48* 


3.48* 


3.48 


-0.05 


-0.05 


0.05±0.05 


48Ca 


7.26 


6.27 


8.57 


3.52 


3.53 


3.47 


0.23 


0.24 


0.2±0.05 


90Zr 


7.79 


6.99 


8.66 


4.27 


4.28 


4.27 


0.11 


0.12 


0.05±0.1 


208p]2 


7.25 


6.52 


7.86 


5.49 


5.49 


5.50 


0.26 


0.27 


0.16±0.05 



discrepancies between the results and data, the present model provides reasonable results. 
In particular, as in QMC-I, it reproduces the rms charge radii for medium and heavy nuclei 
quite well. 



§4. Summary 



I have reviewed the quark-meson coupling (QMC) model to include quark degrees of free- 
dom in the hadrons involved, and have shown the density dependence of hadron masses in 
^'^Ca as an example. As several authors have suggested [l|, Q, the hadron mass is reduced 
because of the scalar mean-field in medium. In the present model the hadron mass can be 
related to the number of non-strange quarks and the strength of the scalar mean-field (see 
Eq.([Tl|)). We have found a new, simple formula to describe the hadron masses in the medium, 
and this led to a new scaling relationship among them (see Eq.(12)). We should note that 



the origins of the meson- mass reduction in QMC and QHD are completely different |^, IC]. 
It would be very interesting to compare our results with forthcoming experiments on the 
vector-meson mass |^. 

At low energy an effective field theory to describe nuclei will generally contain an infinite 
number of interaction terms, which incorporate the compositeness of the low-energy degrees 



of freedom, namely the hadrons |16], and it is then expected to involve numerous couplings 
which may be nonrenormalizable. Thus, one needs an organizing principle to find out a 
sensible field theory. Manohar and Georgi [|l^] have proposed a systematic way to manage 
such complicated, effective field theories called "naive dimensional analysis" (NDA). NDA 
gives rules for assigning a coefficient of the appropriate size to any interaction term in an 
effective lagrangian. After extracting the dimensional factors and some appropriate counting 
factors using NDA, the remaining dimensionless coefficients are all assumed to be of order 
unity. This is the so-called naturalness assumption. 

Our lagrangian density, Eq. dH) , provides a lot of effective coupling terms among the meson 
fields because the mesons have structure. In particular, the lagrangian automatically offers 
self-coupling terms (or non-linear terms) with respect to the a field. Therefore, it is very 



interesting whether the QMC model gives natural coupling constants. In Ref. |15], we exam- 



ined the QMC model using NDA, and concluded that the QMC model is quite natural as an 
effective field theory for nuclei. 



The author thanks Tony Thomas and Kazuo Tsushima for valuable discussions and com- 
ments. 
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